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Abstract 

A rotor-router walk is a deterministic version of a random walk, in which the 
walker is routed to each of the neighbouring vertices in some fixed cyclic order. We 
consider here directed covers of graphs (called also periodic trees) and we study sev- 
eral quantities related to rotor-router walks on directed covers. The quantities under 
consideration are: order of the rotor-router group, order of the root element in the 
rotor-router group and the connection with random walks. 
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1 Introduction 

Given a finite connected and directed graph G, one can construct a labelled rooted tree 
T in the following way. The root vertex is labelled with some i S G. Recursively if a; is a 
vertex in T with label i € G, then x has dij successors with label j. The tree T is called 
the directed cover of G. Random walks on directed covers of graphs have been studied 
by Takacs }Tak97j . Nagnibeda and Woess |NW02| . On infinite graphs, their methods 
have been extended by Gilch and Muller [GMllj. 

Rotor-router walks are deterministic analogues to random walks, which have been first 
introduced into the physics literature under the name Eulerian walks by Priezzhev, 
D.Dhar et al [PDDK96J as a model of self organized criticality. In a rotor-router walk 
on a graph, equip each vertex with an arrow (the rotor) pointing to one of the neighbours 
of the vertex. A particle performing a rotor-router walk carries out the following proce- 
dure at each step. First it changes the rotor at its current position to point to the next 
neighbour, in a fixed order chosen at the beginning, and then moves to the neighbour 
the rotor is now pointing at. These deterministic walks have gained increased interest 
in the last years, and in many settings there is remarkable agreement between the be- 
haviour of rotor-router walks and the expected behaviour of random walks. Holroyd 
AND Propp |HP10j proved that many quantities associated to rotor-router walks such as 
normalized hitting frequencies, hitting times and occupation frequencies, are concentrated 
around their expected values for random walks. See also Cooper AND Spencer [CS06J, 
Doerr and Friedrich [DF06j, Angel and Holroyd [AHllj . Kleber [KleQ5], and 
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also Cooper, Doerr et al. [CDFS06J. On the other hand, rotor-router walks and 
random walks can also have striking differences. In questions concerning recurrence and 
transience of rotor-router walks on homogeneous trees, this has been proven by Landau 
and Levine |LL09j . For random initial configurations on homogeneous trees, see Angel 
AND Holroyd [AH IT]. In our note [HS12] we have extended their result to rotor-router 
walks on directed covers of finite graphs. Furthermore, one can use rotor-router walks 
for solving questions regarding the behaviour of random walks: for instance, in [HSllJ wc 
have used a special rotor-router process in order to determine the harmonic measure, that 
is, the exit distribution of a random walk from a finite subset of a graph. 

In this work, we continue the study of several quantities related to rotor-router walks 
such as the order of the rotor-router group and the order of the root element in the 
rotor group on directed covers of finite graphs. On homogeneous trees, this was done by 
LEVlNE [Lev09| . The remainder of the paper is structured as follows. In Section [2] we 
briefly review the definitions and basic properties of: graphs and trees, directed covers of 
graphs, rotor-router walks, rotor-router and sandpile groups and the connection between 
them. We will follow the notation from [HS12]. 

Section [3] is dedicated to the study of the rotor-router group on finite pieces of directed 
covers T of finite graphs G. In particular we study the rotor-router group on balls with 
respect to the graph metric, that is on T h = {x G T : d(r,x) < h}, where r is the root 
vertex of the graph and d(r, x) is the length of the shortest path from r to x. To the graph 
we add a global sink vertex, which is connected to the root and all leaves. By counting 
a certain family of rooted spanning forests, we give a recurrence formula in Theorem 13.31 
for the order of the rotor-router group on T h in terms of the respective orders on the 
principal subbranches of T . Furthermore, in Theorem 13.71 we show that the order grows 
doubly exponential in h. The growth depends on the spectral radius of the adjacency 
matrix of G. Then, we consider the order of the root element in the rotor-router group, 
which can be defined as the number of particles needed at the origin of T h , such that after 
performing a rotor-router walk and stopping the particles when they hit the sink, we are 
back to the same rotor-router configuration we started from. We describe in Theorem 13.91 
a recursive way for finding the order of the root element in terms of the respective orders 
of the subbranches. A key tool in the proof of this result is the so-called explosion formula 
introduced by Angel and Holroyd [AH11| . 

2 Preliminaries 

Graphs and Trees. Let G = (V, E) be a locally finite and connected directed multi- 
graph, with vertex set V and edge set E. For sake of simplicity, we identify the graph G 
with its vertex set V, i.e., i £ G means i 6 V. It will be clear from the context whether 
we are speaking about vertices or edges. If is an edge of G, we write i ~g j- We 
write d(i,j) for the graph distance, that is, the length of the shortest path between i and 
j. A directed graph G is strongly connected if for every two vertices i and j there exists a 
directed path from i to j and a directed path from j to i. 

Let D = (djj)jj e G be the adjacency matrix of G, where dij is the number of directed edges 
connecting i to j. We write di for the sum of the entries in the i-th row of D, that is 
di = J2jeG dij is the degree of the vertex i (or the number of the outgoing edges from i). 
The adjacency matrix D is irreducible if for every pair of indices i and j there exists a 
natural number n such that dff > 0, where d^' represents the (i,j)-entry of the matrix 
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power D n . If G is strongly connected, its adjacency matrix is irreducible. 

A tree T is a connected, cycle-free graph. A rooted tree is a tree with a distinguished 
vertex r, called the root. For a vertex x £ 7", denote by |x| the height of x, that is the 
graph distance from the root to x. For any positive integer h, define the truncated tree 
T h = {x e T : \x\ < h} to be the subgraph of T induced by the vertices at height smaller 
or equal to h. 

For a vertex x £ T \ {r}, denote by x^ its ancestor, that is the unique neighbour of x 
closer to the root r. It will be convenient to attach an additional vertex n°) to the root 
r, which will be considered in the following as a sink vertex. Additionally we fix a planar 
embedding of T and enumerate the neighbours of a vertex s 6 Tin counter-clockwise order 
(x^°\x^\ . . . , x^ -1 )) beginning with the ancestor. We will call a vertex y a descendant 
of x, if x lies on the unique shortest path from y to the root r. A descendant of x, which is 
also a neighbour of x, will be called a child. A cone C x rooted a x is the subtree spanned 
by the descendants of x. The principal branches of T are the cones rooted at the children 
of the root r. 

The wired tree T h of height h is the multigraph obtained from T h by collapsing all leaves, 
i.e. all vertices y £ T with d(r, y) = h, together with the ancestor r^ of the root to a 
single vertex s, the sink. We do not collapse multiple edges. 

Directed Covers of Graphs. Let now G be a finite, directed and strongly connected 
multigraph with adjacency matrix D = (dij). Let m be the cardinality of the vertices of 
G, and label the vertices of G by {1, 2, . . . , m}. 

The directed cover T of G is defined recursively as a rooted tree T whose vertices are 
labelled by the vertex set {1,2,..., m} of G. The root r of T is labelled with some i £ G. 
Recursively, if x is a vertex in T with label % £ G, then x has dij descendants with label 
j. We define the label function r : T — > G as the map that associates to each vertex in T 
its label in G. The label t(x) of a vertex x will be also called the type of x. For a vertex 
x £ T, we will not only need its type, but also the types of its children. In order to keep 
track of the type of a vertex and the types of its children we introduce the generation 
function \ = (x«)jgG with Xi '■ {1> • • • >^«} ~~ * G. For a vertex x of type i, Xi(k) represents 
the type of the fe-th child x^ of x, i.e., 

if t(x) = i then = t(x^), for A; = 1, ... , <ij. 

As the neighbours (x(°), . . . ,x (d -(-) ) ) of any vertex x are drawn in clockwise order, the 
generation function \ & l so fixes the planar embedding of the tree and thus defines T 
uniquely as a planted plane tree. 

In order to distinguish between the two graphs G and T we use the variables i,j for vertices 
in G (and labels or types in T) and x, y for vertices of T ■ The tree T constructed in this 
way is called the directed cover of G. Such trees are also known as periodic trees, see 
Lyons [LP], or trees with finitely many cone types in Nagnibeda and Woess [NW02j . 
We say that the cone C x has cone type r(x). Note that if x,y £ T have the same label, 
that is t(x) = r(y), then the trees C x and C y are isomorphic as rooted trees. Since G is a 
finite graph, the number of isomorphism classes of C x , x £ T is finite. 

The graph G is called the base graph or the generating graph for the tree T. We write 
71 for a tree with root r of type i, that is r(r) = i. In the following we give two basic 
examples of directed covers of graphs. 
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Example 2.1 (Fibonacci tree). The Fibonacci tree is the directed cover of the graph 
G on two vertices {1,2}, with adjacency matrix 




It is a tree with two cone types: a vertex with label 1 (of type 1) in the tree has only one 
child with label 2 and a vertex of type 2 has one child of type 1 and one child of type 2. 
In Figure [1] we have a wired Fibonacci tree with root type 2 (the vertices of type 1 are 
coloured in blue and those of type 2 in red). The generation function x 1S also given in 
the picture above. 

Example 2.2 (Bi-regular tree T). The bi-regular tree T with parameters a,/3 £ N is 
the directed cover of the graph G on two vertices {1,2}, with adjacency matrix 

*-G o) 

It is a tree with two cone types: every vertex in T with label 1 has no child with label 
1 and a = dyi children with label 2, and every vertex with label 2 has j3 = di\ children 
with label 1 and no child with label 2. Since in this case, on each level there are vertices 
of only one type, the function x has to be: Xl(&) = 2, for k = 1, . . . , a and X2(k) = 1 for 
k = l,...,p. 
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2.1 Rotor-Router Walks 

On a locally finite and connected graph G, a rotor-router walk is defined as follows. For each 
vertex x G G fix a cyclic ordering c(x) of its neighbours: c(x) = (x^ ', x^\ . . . , x^ dx ~ 1 ^ , 
where x ~g x^' for all i = 0, 1, . . . , d x — 1 and 4 is the degree of x. The ordering c{x) 
is called the rotor sequence of x. A rotor configuration is a function p : G — >• G, with 
~g for all x G G. Hence p assigns to every vertex one of its neighbours. By 
abuse of notation, we write p(x) = i if the rotor at x points to the neighbour x®, with 
i € {0,1,..., 4-1}- 

A rotor-router walk is defined by the following rule. Let x be the current position of the 
particle, and p(x) = i the state of the rotor at x. In one step of the walk two things 
happen. First the position of the rotor at x is incremented to point to the next neighbour 
aj(*+l) in the ordering c(x), that is, p(x) is set to i + 1 (with addition performed modulo 
d x ). Then the particle moves to position x^ 1+1 \ The rotor-router walk is obtained by 
repeatedly applying this rule. 

Suppose now that G is a finite graph with m vertices and fix a vertex s in G, which will 
represent the sink. 

Rotor- Router Group. Given a rotor configuration p on G, write e x {p) for the rotor 
configuration resulting from starting a particle at x and letting it perform a rotor-router 
walk until it reaches the sink s. If a particle visits a vertex infinitely often, it also visits all 
of its neighbours infinitely often; since G is connected and finite, the particle eventually 
reaches the sink. 

The set of edges {(x,p(x)) : x G G\{s}} in a rotor configuration forms a spanning 
subgraph of G in which every vertex except the sink s has out-degree one. If this subgraph 
contains no cycles, we call it an oriented spanning tree of G. Write Rec(G) for the set of 
oriented spanning trees of G, which is also called the set of recurrent configurations. It is 
easy to see that if p G Rec(G), then also e x (p) G Rec(G). For a proof, see Landau and 
Levine [LL091 Lemma 2.1]. Another interesting property of the rotor configurations is 
that if p\,p2 G Rec(G) and e x {p\) = e x (p2), then p\ = p2] see once again |LL09] for a 
proof. This means that the operation e x of adding a particle at x and routing it to the 
sink acts invertibly on the set of recurrent rotor configurations. 

The rotor-router group RR(G) of G is defined as the subgroup of the permutation group of 
Rec(G) generated by \e x : x G G \ {s}}. For any two vertices x and y, the operators e x 
and e y commute. This is the so-called abelian property of rotor-router walks. Hence the 
group RR(G) is abelian. Furthermore, RR(G) acts transitively on Rec(G). More details can 
be found in [LL09] . 

Sandpile Group. A chip configuration a on G, also called a sandpile on G, is a vector in 
Z m_1 of non-negative integers indexed by the non-sink vertices of G, where o~(x) represents 
the number of chips at the vertex x. A chip configuration a is called stable if o~(x) < d x , 
for every non-sink vertex x. A vertex x is unstable if o~(x) > d x . An unstable vertex 
may topple, by sending one chip to each neighbour. If a is not stable then one can show 
that by successively toppling unstable vertices, in finitely many steps we arrive at a stable 
configuration a° . A stable chip configuration a is called recurrent if there exists a nonzero 
chip configuration 5 such that (<r + 5)° = a. The sandpile group SP(G) may be thought of 
as the set of recurrent chip configurations under the operation (a + 5)° of addition followed 
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by stabilization. The order of the sandpile group SP(G) is given by the determinant of the 



The reduced Laplacian is obtained by deleting from the Laplacian matrix A of G the row 
and the column corresponding to the sink. By the matrix-tree theorem, this determinant 
equals the number of oriented spanning trees of G rooted at the sink. 

Theorem 2.3. The rotor-router group RR(G) for a connected finite graph G with a global 
sink is isomorphic to its sandpile group SP(G). 

The proof can be found in [LL091 Theorem 2.5]. 

Notation: two non-negative functions f{n) and g(n) have the same growth and we write 
f(n) >c g(n) if there exist constants c\,C2 such that c\g{n) < f(n) < C2g(n). 

For the rest, we fix the following: 

• G finite graph with m vertices labelled by {1,2,... , m}. 

• Ti directed cover of G with root r of type i G G. 

• 7^ the wired directed cover of height h and root type i. 

• t(x) label function and Xi(k) generation function. 

3 Order of the Rotor-Router Group 

In this section we want to describe the rotor-router group on directed covers of finite 
graphs. For homogeneous trees, this was done in LEViNE jLev09j . and the method used 
there fails when one considers non- homogeneous structures. Because of the non homo- 
geneity, our approach is also quite technical. We will relate the rotor-router group of a 
wired directed cover of a graph with the rotor-router group of its principal subbranches. 
Denote by t% the number of spanning trees of T/ 1 . 

From Section \2.1\ in order to find an expression for the order |RR(7^' l )| of the rotor-router 
group of 77 1 , it is enough to count the number of the spanning trees, since | RR(l~ h )\ = t\. 
We will count spanning trees in terms of a class of spanning forests of the original tree 
l~ h , which we define next. 

Definition 3.1. We say that a spanning forest of a graph H is rooted at a set of vertices 
S C H if every connected component of the forest contains exactly one vertex of S. 

Denote by H/S" the multigraph obtained from H by contracting S into a single vertex s, 
while keeping multiple edges that may have been created by this process. 

Lemma 3.2. There is a bijection between the set of spanning trees o/H/5 and the set of 
spanning forests of H which are rooted at S. 

Proof. Let F be a spanning forest of H and denote by F j S the spanning subgraph of H/S 1 
obtained by contracting S into a single vertex s. If F has a connected component which 
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does not contain a vertex of S, its contraction F/S is still not connected. On the other 
hand, if F has a connected component which contains at least two vertices of S, collapsing 
S into a single vertex creates a cycle in F/S. Hence F/S is a spanning tree of H/S if and 
only if F is a spanning forest rooted at S. 

Let Eg be the set of edges e of H such that not both endpoints of e are contained in S. 
Consider the map tp : H — > H/S, defined by 



ip{v) 



s, for v E S 
v, otherwise 



and its natural extension to the edge set E of H which we also call tp. Then ip is a bijection 
between Es and the edges of H/S. Since every spanning forest F of H rooted at S is a 
subset of Es, the map ip naturally extends to a map t/j : F \— > F/S, which maps F onto 
its contraction. Hence tp is injective, since the tree F/S is fully defined by its edge set. 

Now let T be a spanning tree of H/S, and define F = {tp~ 1 (e) : e is a edge of T}. Then F 
is a forest from which we may obtain a spanning forest rooted at S by adding every vertex 
of S which is not contained in a connected component of F as a single vertex component. 
This construction is again injective. Hence there exists a bijection between the set of 
spanning trees of H/S and the set of spanning forests of H rooted on S. □ 

Back to the truncated tree F/ 1 , let us introduce the down and up sinks 

Sl = {A ^} and s r = {x G T/ 1 : \x\ = h}, (1) 

where n°) is the ancestor of the root r, an additional vertex connected with the root r of 
T/ 1 , and let 

S = s± U s r . (2) 

With the notation introduced above we have T/ 1 = F/ 1 / S. According to Lemma 13.21 in 
order to compute the order of the rotor-router group, one has to count the number of 
spanning forests of T/ 1 rooted at S. We partition the spanning forests into two types. For 
all i £ G denote by 

Fj 1 , = the number of spanning forests of 7^ rooted at S, 
which contain an edge from r to s±, 

F^ = the number of spanning forests of F/ 1 rooted at S, 
which contain a path from r to s . 

Using this notation we have that 

\RR(jf)\=t^ = F^ + F^. (3) 

and the order of the rotor-router group can be calculated recursively as follows. 

Theorem 3.3. Fhe order of the rotor-router group RR(7^ ft ) is given by |RR(7^ )| = 

F-, + FK, and the number of spanning forests Fj 1 ^ and FK with i G G can be calcu- 
lated recursively as: 

d .. F h-i (4) 

j£G r j,l +jt j,t 

The initial values are Fh = 1 and FK = d{, for all i G G. 
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Proof. Let T/ 1 be a tree with root r of type i € G, and of height h. For k = 1, . . . , di, we 
denote by T^7^\ the cone of the fc-th child r( fc ) of the root r in 7^. By the construction of 

the directed cover, the root of 7^.7^ is of type Xi{k)- 

Every spanning forest J- that is counted in F h , contains the edge (r, sj,). In this case, 

2,4- 

whenever the edge (A k \r) is contained in T , its restriction to T h 7u\ does not contain a 

path from to the upper sink of T^T^y On the other hand, if the forest T does not 

contain the edge (A k \r), then the restriction of the spanning forest to the first level cone 
must contain a path from A k ^ to the upper sink. Otherwise the connected component 
containing r^) would not contain a sink vertex. As we can freely choose which of the 
vertices {(A k \ r) : k = 1, . . . , di} are part of the forest, we get 

di 

where w = (u>i,u>2, . . . , w^) is a word of length di over the alphabet {4-,t}. In order to 
reduce the previous equation to the form given in (£[]), take a j 6 G with dij ^ 0. This 
implies that for some k = 1, . . . , di we have Xi(k) = j (this happens exactly dij times) and 
FJ 1 , can be factorized as 

weRt}^- 1 fe=i 

This procedure can be repeated exactly <ij 7 - times, and we obtain 



F h , = [F h 7 1 + F^ 1 ) ° Y TT F h -\ 
By proceeding in the same way for all j € G for which dij ^ 0, we get 

p t= n i^ 1 + F ^T j , 

which can be obviously extended to all j 6 G, and we have proved the first part of @). 

Consider now a spanning forest T with the root connected to the upper sink s . Thus, 
whenever m edges from the set {(A k \r) : k = 1, . . . ,di} are contained in J 7 , there are 
exactly di — m possible ways to connect r to s . Hence 

*?t= e ^wkj,, (5) 

w e{i,\} d i k=i 

where U(w) = #{k G {1, . . . , di} : Wk = t} represents the number of letters t in the word 
u;. Note that the number of spanning forests depends only on the adjacency matrix D of 
the generating graph, and not on the particular planar embedding \ of the tree. Therefore, 
to prove the equivalence of ^ and 0, we can use induction on the number of children 
of a vertex of a certain type. 
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For the rest of the proof, fix i G G and denote by D % the i-th row of the adjacency matrix 
D of G. If a vertex of type i has only one child, say of type I, that is, D l = (^'(0)j 6 G' 
then ([5]) reduces to 

*ft= £ U{w)F^ = F^ 
«>e{4.,t} 

which proves the induction base. Assume now that for all possible choices of F.^ 1 , F^ 1 G 
M and of the adjacency matrix D, the following holds 

Me {4.,t} d i fc=i jeG jeGj,4 " r j,t 

and we prove that it also holds if we increase one entry of D by 1. Let us suppose that 
we increase the dj;-entry by 1, for some fixed I G G. This means that any vertex of type i 
has an additional child of type I. Let us denote by D this new matrix and by d, and Xi 
the quantities corresponding to the degree of a vertex and to the generation function for 
D. Due to the fact that the number of forests does not depend on the planar embedding, 
i.e., on the function %j, we can put the new additional child of type I on the rightmost 
position, that is Xi(di) = I and Xi(k) = Xi(k) for k < dj. Then we have 

e %)k; t = E u ^(uka^ f ^i +1 

we{l,1-} d i k=l w£{l,t} d i +1 \fc=i / 

which, factoring the word w by its last letter, is equal to 

= E w«>+i)(n j SiuW I + E ^wflk^wW 1 
= (^r 1 +^r 1 ) E ^n^+^IK^+^f 1 )^- 

we{4-,t} d i fc=i jeG 
Using the induction hypothesis (JSJ) this further equals 

n (itr 1 + pjr 1 )^® E + F tr l n + ^f 1 )^ 

jeG jeG j,4 + j,t jeG 

nfph-l , rA-l\i«+i((i) \ - (jg +^(0) i tt 1 
^ j,4- + j,t J ph-l , F h-1 ' 

jeG jeG j,4- " l j,t 

which proves the inductive step. □ 

Next, we evaluate asymptotically the behaviour of |RR(7^ ?i )| = F^ + F^, for large values 
ofh. 

Lemma 3.4. For i G G and li£N, £/ie sequence y£ defined by 

h _ *i,t 
i,4 

25 coTivcvgcfit. If the spcctrul radius t^D^j > 1 then the limit = lini/j—^oo 7^ is positive. 
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Recall that the spectral radius f(A) of a square matrix A is the maximal absolute value of 
all the eigenvalues of A. To prove Lemma 13.41 we will need two additional results. 

Proposition 3.5. Let A G jj mxm fr e a non-negative matrix and let v G R m be a positive 
vector. If there is a number b > such that Av < bv, then r{A) < b. 

For the proof see Ding and Zhou [DZ09, Proposition 2.2]. The next result is a sim- 
ple application of Tarski's Fixed point Theorem for increasing functions on lattices, see 
Kennan [KenOll Theorem 3.3]. 

Theorem 3.6. Suppose f is an increasing and strictly concave function from R m to R m 
such that /(0) > 0, f(a) > a for some positive vector a, and /(b) < b for some vector 
b > a. Then f has a unique positive fixed point. 



Proof of Lemma \3.4\ Let us fix i G G for the rest of the proof. The initial value is = di 
Substituting in the definition of 7^ the recurrence relations (dJ) for and F- 1 ^, we get 



7 h = y 1 oA = y d h for h e N- (7) 

h F h-i + F h-i Z_, v x + h-i v j 

In order to prove that 7^ is convergent in h, we show that 7^ is monotone and bounded. 
We first show by induction on h that 7^ is strictly decreasing. We know that 7* = di, and 
the induction basis 7? < 7* follows from 

il = Yl ^r+V < 2 dij = di = 7i1- 

jeG ^ h jeG 

For the inductive step, we assume that 7^ < %~ , and we prove that 7^ +1 < 7^. This 
follows from 

The sequence 7^ is also bounded, i.e., < 7^ < di, therefore it converges to a limit. 
Denote by Yi = lim^oo 7i an d T = (T\, . . . ,T m ) G K> - Then the limit vector T is a 
solution of the fixed point equation 

^ = E^rTr- ( 8 ) 

jeG ^ ' 

In order to prove the positivity of the limit vector T we apply Theorem 13.61 to the function 
/=(/i,.. .,f m ) :R m ^M. m with 

fi(x) = ^dij Xj and x = (x 1} . . .,x m ). 

*-t 1 + Xj 

3 eG J 

The function / is obviously strictly concave and increasing. From the first part of the 
proof we have /(b) < b for b being the vector of initial values b = (d\, . . . , d m ). 

Let now r(D) > 1 and assume that for all positive vectors a = [a\, . . . ,a m ) we have 
/(a) < a, which can be written as 



/ a 

/ l+o: 



< a. (9) 



l+a„ 
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Let a* = maxjaj : i = 1, . . . , m}, then © implies ^j—* Da < a, hence Da < (l+a*)a, and 
from Proposition 13, 51 it follows that r(D) < 1 + a*. Since a was arbitrary, a* can be made 
arbitrary small, thus r(D) < 1 which is a contradiction to our assumption. Therefore there 
exists a positive vector a = (a±, . . . ,a m ) such that f(a) > a. Theorem 13.61 now ensures 
the existence of a unique positive fixed point of /, hence the equation © has a unique 
positive solution in addition to the trivial solution. 

It remains to show that the limit vector cannot be zero. Let 7^ be the sequence defined 
as 7^ in ([7]) with initial values 7* = a^. The vector a is such that f(a) > a. It is easy 
to see that the sequence 7^ is increasing and 7* — 7 1 > 0. By induction on h, supposing 



li ~ li > 0, we get 



h+l 



h+l 



1 + ^ 



1 + 7* 



5>. 



(l+7f)(l + 7f) 



>0, 



therefore 7^ — 7^ is non-negative. Hence 7^ is bounded from below by 7^, and the limit 
Ti is positive. □ 

Theorem 3.7. Let G be a finite, directed and strongly connected graph with vertex set 
{1,2,... , to} and adjacency matrix D. Let % be the directed cover of G with root r of type 
i and T* be the wired tree defined as above. Ifr(D) > 1, then the order of the rotor-router 
group RR(7^ ) grows doubly exponential: 



|RR(7f)| x exp{r(D) h }, for all i € G, 
where r(D) is the spectral radius (Perron-Frobenius eigenvalue) of D. 



(10) 



Proof. Recall first that |RR(7^' l )| = F^ + F^. In order to simplify the system of equations 
dH) let us make the following substitution: for i £ G denote by x^ = log and = 
log F^ and apply the logarithm function to (|3J . We get 

^=E^r i +E^ i °g( i +^ 1 ) 

jeG jeG 
yf = 4 + log ( 7 f). 

Consider now the following vectors in M. m : 



x 



fx h \ 


, y h = 


(y h \ 


\, i h = 








\y h m J 







Then (1111) can be written in matrix form as 



x h = Dfx*- 1 +log(l + 7^)) 
y fc = a J ft + log (7"), 



fill 



(12) 



where the function log in (|12p is applied componentwise to the entries of the vectors ~y h and 
(1 + 7' 1 ) respectively. The initial values are x 1 = (0, . . . ,0) T , y 1 = (logdi, . . . ,logd m ) T , 
and 7 1 = (di, . . . , d m ) T . Then the solution of (fT2|) is given by 



h-l 



x h = Y J D k log(l + 7 7i - fe ) 



(13) 



k=l 
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From Lemma 13.41 the entries of the vector log (l + -y h fc ) are bounded: for all i G 
{1, 2, . . . , m} and k = 1, . . . , h 

< log(l + Ti) < log(l + jt") < log(l + 7i) = log(l + di). 
All Ti and di are positive. Write q = log(l + Tj) and Cj = log(l + dj). Then 

(ij^)^ 1 *^ / c ' (14) 

where c = (ci, . . . , c m ) T and C = (C%, . . . , C m ) T . The sum 

(Eti behaves like D h 

for big values of /i, and the exponential growth rate of the matrix power D h as h —> oo 
is controlled by the eigenvalue of D with the largest absolute value. Since D is a non- 
negative and irreducible matrix, according to Perron-Frobenius theorem for irreducible 
matrices, there exists a positive real number r{D) (the spectral radius of D), called the 
Perron-Frobenius eigenvalue which is the eigenvalue of D with the largest absolute value. 
By (|14p . we can conclude that the behaviour of x h , for h — > oo is given by the greatest 
eigenvalue r(D) > 1 of D, i.e., for all i 

4 x r{Df. 

Since y h differs from x h only by a bounded and decreasing quantity, see (|12p. there exists 
c > such that x r(D)' 1 + c. Using = log Ft 1 , and = log F^, we get 

ift x expMD)' 1 } and ift x exp{r(D) h }, 

which implies 

|RR(^)|xexp{r( J D)' 1 }, 
and this proves the statement. □ 

Example 3.8 (The Fibonacci tree). For the Fibonacci tree, the system of equations §S§ 
can be written as 



rph _ jph-1 , ph-l 

jph _ jph-l 
r l,t ~~ 2,t 



F 2 h f = F^ (F^ + F^) + F^ (F^ + F^ > 



i,t 



The two sequences 7^ and 72 are given recursively by 

„,>i-l ../i-l 
7i = - , h _i and 7 2 = - + 



The initial values are 7* = 1, 7^ = 2 and the limit values are 

Ti = lim ^ = V2-1 and T 2 = lim 7^ = — . 

h— >oo h— >oo 2 

The Perron-Frobenius eigenvalue of D is , which is also related with the Fibonacci 
numbers F n by limn-^ F g +1 = l+ y^ . Finally, regarding the order of the rotor-router 
group RR(7^' 1 ) we have 



\RR(T?) 



cxp 
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3.1 Order of the Root Element in the Rotor-Router Group 

Levine |Lev09| computed the order of the root element in the rotor-router group on 
homogeneous trees. His approach holds only for homogeneous trees, and it fails in our 
case. We describe here the order of the root element in terms of the respective orders on 
the principal subbranches. Our method is yet another useful application of the explosion 
formula introduced in [AHlll Theorem 11]. 

Like before, 71 is the directed cover of G with root r of type r(r) = i. The principal 
branches of % are the subtrees 7~ Xi (k) rooted at the children A k ^ of the root r with type 
Xi{k) = r(A k )) and k = 1,2,... Hence % has as principal branches di = SjeG ^ij 
subtrees. Recall the definition of the truncated wired tree 7^ , with the same root r and 
sink S = S| U s^. Suppose now that we have one particle at the root r of T^ 1 , and we 
let it perform a rotor-router walk until it hits the sink S, where it stops. Denote by fh 
the element of the rotor-router group RR(7; ) corresponding to this process. Similarly, 
for all k = 1 , . . . , di , denote by the element of the rotor-router group RR(7^7 fe \) 
corresponding to one particle performing rotor-router walk on the principal subbranch 
starting at A k \ Write (fh) for the cyclic subgroup of RR(7^ /l ) generated by fh, and 

{{^h-v • • • >^i-i)) fo r the cyclic subgroup of ©f l =1 ^C^wifc 1 )) generated by the element 

(fh—ii ■ ■ ■ ' ^i— l) • ^ |Lev09l Theorem 3.3] and from the isomorphism between the sandpile 
and the rotor-router group of a tree, we have 

k=l 

For simplicity of notation, we denote by R% the order of the element fh in the rotor-router 
group RR(7^' 1 ), that is, the cardinality of the cyclic group (f h ), i.e. R!> = \ (f h )\. Recall the 
definition of the down and up sinks sj, and s^ respectively and let 

= the number of particles stopped in s^ after R^ particles 
have been routed from r, 

= the number of particles stopped in s t after R% particles 
have been routed from r, 

for i € G. Then 

Rt = S^ + S^. (16) 

Write R ~}\ for the number of particles started at the origin A k ^ of the branch T h 7,\- 

Theorem 3.9. The order R% of the root element in the rotor-router group is given by 
RJl = + and S^, can be computed recursively as follows: 

S i,l = lc ™({ R Mk) : fork = l,...,di}) 
ah _ ah 

iec R j 

for all i £ G. Here 1cm represents the least common multiple. The starting values are 
given by = 1 and S^*. = di, for all i e G. 

The proof of this result uses the explosion formula introduced in [AH11, Theorem 11]. 
Since we will need this formula in the proof, we first adapt it here to our case. 
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3.1.1 Explosion Formula 

The explosion formula, introduced in [AHlll Theorem 11], gives a recursive formula for 
computing the number E n (Ti,p) of particles which escape to infinity when we start n 
rotor-router walks at the root r of a tree %■ The initial configuration of rotors on % is p. 
For the tree % and principal branches T Xi nA, rooted at the children of r, write p Xi (k) 
for the restriction of the rotor configuration p on 7wfc), with k = 1, . . . , dj. 

Let us first introduce some notations, following mainly the notations from [AHll j Theorem 
11]. Let e n = e n (Ti,p) = 1 [particle n escapes to infinity] , where 1 represents the indicator func- 
tion, so that E n = Ylk=i e k- Moreover, let e(%, p) be the escape sequence (ei, e 2 , . . .). Let 
N + = {1, 2, . . .} and for sequences in N N+ we denote addition by (oj, a 2 , ■ • .)+(&i, 6 2 > • • •) := 
(01 + 61, a% + 62, • • •)• Define the shift operator 8 by 

0(ai,a 2 ,...) := (0, ai, a 2 , • • •) 

and the explosion operator X by 

£(01,03,...) :=(l a \0,l a2 ,0,...) 

where l k denotes a string of k Is, or the empty string if k = 0. We define a majorization 
order ^ on sequences (ai, a 2 , . . .) ^ (&i, 6 2 , . . .) if and only Yy{ a j ^ Yli bj for all n. For 
a = (ai, a 2 , . . .) denote by a^ = its k-th element. For a finite sequence a the length is 
denoted by |a| and a n is its n times repetition. We now adapt jAHlll Theorem 11] for 
directed covers of finite graphs. 

Theorem 3.10 (Explosion formula). Let G be a finite graph with m vertices and % its 
directed cover with root r of type i. Fix a rotor configuration p on %■ Then 



\ fe=l 



fc=p(r)+l 



Remark 3.11. In case of a finite tree l~ h , we define the escape sequence e such that 
e n (7l h ,p) = 1 if the n-th particle reaches the upper sink s^ before reaching s^, and 
otherwise. Then the explosion formula can be also written in the form 

e(T i \p) = x( P jr9e(T x %),p Xi{k) )+ £ <T^ y p xm )\ 

\ k=l k=p(r)+l / 

with initial values e(7^°, p) = (1, 1, . . .) for all i G G. 

Proof of Theorem \3. 91 Fix an i £ G. It is enough to consider only the zero rotor config- 
uration 0, that is, at all vertices the rotors point to the ancestors. Recall that R% is the 
order of the root element in the rotor-router group of '■ Hence, after R^ particles have 
performed a rotor-router walk (stopped at the sink S) we are back to the configuration 
0. In other words, the escape sequence e(7^\0) is periodic with period Rf. Write e(7^) 
for the finite sequence consisting of the first full period of e(7^,0). In order to adapt the 
explosion formula to the finite sequence e{T^), all the escape sequences involved need to 
be extended such that they have the same length. Let 

lcm (V ( R h () 



p/i 
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We can now write the explosion formula for the first particles as 

e(T^ = x(,r) with vr=£e(7^)^ 1( *>. (18) 

k=l 

Hence eiT/ 1 ) is a sequence of length = + consisting of "zeros" and "ones", with 
being the number of O's and S-* being the number of l's. Since the number of O's in 
a string X(a) is equal to |a|, we get 

^ = |vn=lcm(< {1)) ...,^ } ). 

On the other hand, because equals the number of l's in the escape sequence e(1~ h ), 
the explosion formula (JTBj) gives 

Ri ^ s ti ( ck \ 

1=1 1=1 1=1 \k=l ) i 

In the last equality we used the fact that for any finite sequence a, Yl\=i a ' = Si=i ^( a )i 
holds. By exchanging the order of summation in (|19p . we obtain 

ah R h ~ 1 

k=l 1=1 v 7 1 fc=l i=l 

di ah-1 qh-1 

fc=i fc=i xi(*) jeG 

which proves the theorem. □ 

We give now an alternative way of writing the system (|17p . Using (|16|) , let us add the two 
equations in (fT7|) . and then divide the result through R^. We get 

1 

Because of the 1cm involved in ()17p . it is hard to derive asymptotics for Rf. 

Connection with the Random Walk. Let (Xt) be a simple random walk on 71 which 
starts at the root r of H and define the stopping times Tj (which depend on the root r 
with type r(r) = i) as 

T< = inf{t >0:Iie Si u/} 
and hitting probabilities ± and if j f as 

H iA = P r [X Tt = Sl ] and = ¥ r [X Ti € s f ] (21) 

Write U?\ and iJ. ^. for the corresponding probabilities on the truncated tree T- 1 . It is easy 
to see that H it ± (or if^t) can be expressed as the quotient between the number of particles 



lcm <^ R 



({■ 



? h-l 



for 



i?' 



di + 1 - d 



S 



i,4- 



(20) 
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l ( or ^if) that are routed in the sink s± (or s ) and the total number of particles R\ 
started at the root, that is 



H, 



oh 

ft? 



and H 



t 



oh 

°i,t 



One can also obtain that Hi is the solution of the same equation (|20j) . by factorizing the 
random walk (Xf) with respect to the first step 



d, 



i = hi (di + 1 - Km+) = H l + 1 - E ^r 1 ) ■ 
fc=i j=i 



(22) 



with the initial value iT . , = l/(di + 1). 



Example 3.12 (Fibonacci Tree). Consider the order of the root element in the rotor- 
router group RR(7^' 1 ) of the wired Fibonacci tree T/ 1 . In this case the equation ([2Uj) is 
given as: 



oh 



R. 



h-l 



Rl = S x \ 2 



S: 



h-l 

2,1 



R: 



h-l 



h-l r>h-l\ 



S 2 \ =lcm {Rr\R n 2 



S. 



s 



h-l 

1,1 



S 



ft.-l ' 

2,1 



2.1 



R 



h-l 



R 



h-l 



Computations suggest that on the Fibonacci tree, the sequences S^ and R% are relatively 



prime. 



Example 3.13 (Bi-regular tree). For the (a, /3)-bi-regular tree the recurrence relation 
(|20p for the order of the root element reduces to the simple form 

R^ = R%-\a + l) -aR^ 2 

R% = RtHP + i)-mt 2 , 



which has the explicit solution 

' (a/3) fc - 1 

R\ = 



af3 - 1 
(a/3) fc+1 - 1 



a(3 - 1 



{a + 1)13 + 1, forh = 2k 
a + 1), iovh = 2k+l. 



For i?2 the solution is the same, with the roles of a and f3 exchanged. 
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